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Abstract 

This paper describes the data analysis technique used for magnetic testing at the NASA Goddard Space Flight 
Center (GSFC). Excellent results have been obtained using this technique to convert a spacecraft’s measured magnetic field 
data into its respective magnetic dipole moment model. The model is most accurate with the earth’s geomagnetic field 
cancelled in a spherical region bounded by the measurement magnetometers with a minimum radius large enough to enclose 
the magnetic source. Considerably enhanced spacecraft magnetic testing is offered by using this technique in conjunction 
with a computer-controlled magnetic field measurement system. Such a system, with real-time magnetic field display 
capabilities, has been incorporated into other existing magnetic measurement facilities and is also used at remote locations 
where transport to a magnetics test facility is impractical. 
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Introduction 

A spacecraft whose magnetic field has been accurately calibrated possesses a unique capability for acquiring 
magnetic field and other scientific data. Geomagnetic, as well as planetary and interplanetary magnetic measurements depend 
on such a controlled spacecraft magnetic environment for accurate results. To provide this environment, restraint of magnetic 
materials on board the spacecraft and sophisticated prelaunch magnetic tests are required. Some facilities, like the GSFC 
Spacecraft Magnetic Test Facility (SMTF), are able to null Earth’s geomagnetic field to less than 0.5 nT, providing a 
controlled environment for measurement and control of spacecraft fields. 

The magnetic field modeling technique described herein was developed at the Goddard Magnetic Test Facility as an 
aid in the evaluation of spacecraft magnetic fields and has been used for over thirty years. Three specific problem areas 
requiring magnetic field measurements have been investigated using this technique: 

• Spacecraft attitude control problems resulting from a large intrinsic magnetic dipole moment; 

• Reduced accuracy of data returned from magnetometer experiments due to spacecraft bias fields; and 

• Complicated data analysis caused by spacecraft magnetic fields around low-energy electron experiments. 

The magnetic field modeling technique described in this paper, combined with a computer measurement system, provides a 
means of solving the above problems. Using the computerized modeling technique to provide a near-real-time visual 
representation of the spacecraft magnetic field reduces testing time required for complex magnetic measurements. This visual 
representation also aids in reduction of the dipole moment by showing the correct orientation for placement of compensation 
magnets or demagnetization treatment. In addition, the magnetic field modeling technique provides a rigorous method of 
determining the magnetic field at any point of interest around the spacecraft, based upon a limited number of measurements 
at discrete points. This enables a rapid and comprehensive evaluation of spacecraft magnetic fields for the solution of 
magnetometer bias and low-energy electron deflection problems. 



Background 


In 1972, William L. Eichhom published a NASA technical report describing a dipole moment measurement 
technique using near-field magnetic data [1]. The technique assumes that a spacecraft’s magnetic field can be represented by 
the addition of a number of multi-pole fields. These fields are measured in the near-field region that begins on a spherical 
surface, which totally encloses the magnetic source and extends outward until the fields become dipolar. Near-field analysis 
allows dipole moment data to be taken closer to the spacecraft where magnetic fields are strong enough for reliable 
measurements. The technique has been used with good results at the Goddard Magnetic Test Facility for magnetic testing of 
various spacecraft. 

A typical spacecraft orientation for magnetic testing is shown in Figure 1 (where to = direction around axis of 
rotation). Four tri-axial magnetometers are located along one of the spacecraft axes, and magnetic field data is acquired while 
rotating the spacecraft about the vertical (Z) axis past an array of fixed magnetometer probes. The equations for calculating 
the magnetic fields present at the X-, Y-, and Z-axes of a magnetometer sensor located a distance R from the center of the 
spacecraft (as a function of rotational angle 0) were derived in [1]. With a few minor changes in notation and truncation of 
the infinite series at n = 4, to provide a finite series for modeling purposes, these magnetic field equations become: 

B X (R,0,O) = + l )/i?" +2 k m cos md + B„ m sin m£?)C( 0) 

B y (R,0,O) = YLJLo m / Rn+2 ( A ™ sin md + cos m0 ) p n m ( 0) 

B Z (R,0,O) = r n ,i:j Rn+2 (^ cos md + sin m0)d/dt{p:( 0)) 
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magnetic field component (nT) along X magnetometer axis as a function of R and 6; 

radial distance from center of magnetic mass to magnetometer, which must be greater than 
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vertical angle with respect to spacecraft axis of rotation. 


i 

z 



Figure 1. Spacecraft Orientation for Magnetic Testing 



In 1973, the magnetic field modeling technique was developed by Thomas N. Roy to extend near-field analysis. The 
technique is capable of producing accurate mathematical models from the same near-field analysis data. Near-field analysis 
involves only the cos 0 and sin 0 terms of the above equations for calculation of X- and Y-axis dipole moments, and the 
constant term with respect to 0 for the Z-axis moment. Extension of near-field analysis for magnetic field modeling involves 
as many as eight harmonic terms of the magnetic field equations. Experience in using these equations to model fields around 
several different spacecraft has shown a solution order of n = 4 to be sufficient for accurate modeling over regions of interest. 
The technique can easily be extended to higher ordered harmonics for modeling the extreme near-field regions of strongly 
multi-polar sources. 

This paper describes the magnetic field modeling technique and its application to spacecraft magnetic fields. 
Limitations of this technique, as well as other possible uses, are discussed at the end of the paper. 

Magnetic Field Modeling 

Results obtained from solution of the magnetic field equations for n = 1 to 4 indicate that the magnetic field of a 
spacecraft can be described anywhere on a plane containing the measurement magnetometers by: 

Bj(R,#) = DCjR A) +Y > „JC ] n)cos( n d)/R Bn > +Y' n JS j n)sm( n 0)/R Cnj 

Where: 

j = measurement magnetometer axes X, Y or Z; 

R = radial distance; 

6 = rotational angle with respect to magnetometer axes; 

n = order of Fourier solution usually truncated after eight terms; 

DCj = constant field (with respect to 0) modeling coefficient calculated for each axis; 

Cjn = cosine field modeling coefficient calculated for each axis and n = 1 to 4; 

Sjn = sine field modeling coefficient calculated for each axis and n = 1 to 4; 

Aj = constant field decay power calculated for each axis; 

Bnj = cosine field decay power calculated for each axis and n = 1 to 4; and 

Cnj = sine field decay power calculated for each axis and n = 1 to 4. 

Therefore, to construct a magnetic field model using the above equation, a method is needed to calculate the field modeling 
coefficients and field decay powers. 

The procedure for calculating field modeling coefficients and field decay powers begins with this fourth-order 
Fourier analysis. The analysis produces eight Fourier coefficients for each axis (j) (four for the Z-axis) of every 
magnetometer (i), and their standard deviations. Fourier coefficients typically represent averaged values for five spacecraft 
rotations. 

Four cosine Fourier coefficients, calculated for n = 1 to 4: 

Four sine Fourier coefficients, calculated for n = 1 to 4: 

Constant Fourier coefficient: 

■xwSl.s/M 

Where: 

i = magnetometer number, 1 to 4; 

Bji (Rj,0) = rotational magnetic field; and 

j = magnetometer axes X, Y, or Z. 

Next, dividing them by their respective Fourier coefficients normalizes the standard deviations. The normalized cosine 
standard deviations, calculated for n = 1 to 4 are: 


Cjn^aCj^/Cjn, 



The normalized sine standard deviations, calculated for n = 1 to 4 are: 


<r. Sjn, =aSjt %,/Sjtt, 

The normalized constant (Z-axis) standard deviation is: 

<r nJ DC, =aDC Ji /DC Ji 

At this point, the technique assumes that for each axis (j) and each order of Fourier analysis (n), the Fourier 
coefficients from each magnetometer have the same polarity and decrease with increasing radius (R). This permits the 
calculation of field decay powers for each order of Fourier analysis starting with the constant term. All possible combinations 
between the four magnetometers are averaged. Each combination is weighted using the normalized standard deviations 


calculated above: 


Aj 

Bnj 

Chj = 

constant field decay power; 

cosine field decay power calculated for n = 1 to 4; and 
sine field decay power calculated for n = 1 to 4. 


For example, using DC Fourier coefficients and four magnetometers, where R h R 2 , R 3 , and R, are the actual 
magnetometer radii from the center of rotation, calculating the constant field decay power: 

(Aj\ 2 = log[(DCy, )/(DC j2 )]/log(R, /Rj) 

(Aj \ 3 = log[(X>C )/(DC j2 )]/log(R, /R 3 ) 

(Aj \ A = log[( DC j, )/(DC, a )]/log(R, /R 4 ) 

(Aj\ 3 = \og[(DC j2 )/(DC J3 )]/log(R 2 /R 3 ) 

(Aj) 2 , 4 = \og[(DC J2 )/(DC J4 )]/\og(R 2 /R 4 ) 

( Aj\ 4 = log[(DC j3 )/(DC j4 )]/log(R 3 /R 4 ) 

The averaged field decay power is thus: 

[(A J ) 12 /(a„ DC jX )(<T n DC J2 )] + [(Aj ) 2 3 l(a„DC j2 )(a H DC j3 )] + 

[(Aj) x4 KcT'DCj 3 )(<t„DC J4 )] + [(Aj ) v !{a„DC n )(a„DC J3 )] + 

[(Aj), 4 !{a n DC jX )(<r n DC j4 )1 + [(^ ) 24 l(a„DC j2 ){a n DC j4 )] 

A J~ 1 

[(a„DC j{ )(a„DC j2 ) + (a„DC J2 )(<J„DC j3 ) + (a„ DC n )(a n DC j4 ) + 

(a n DCj, )(a n DC J3 ) + (a n DC JX )(a n DC J4 ) + (cr„DC j2 )(ct„DC j4 )] 

When all the field decay powers have been calculated, the Fourier coefficients are normalized using their respective 
magnetometer radii and field decay powers to produce magnetic field modeling coefficients, These field modeling 
coefficients are a result of a weighted average of each set of Fourier coefficients over all four magnetometers using their 
respective normalized standard deviations to provide the weighting factors: 

C J n = Z U C M R i) Bni faW/Z IVr-Cjn, 

Sjn = Zm(V.(*/)° ,J 

dc j = ZliDCAR,) 41 hDC^/'Z'j/^DCj, 


Where: 


Cjn 

Sjn 

DCj 


cosine field modeling coefficient calculated for n = 1 to 4; 
sine field modeling coefficient calculated for n = 1 to 4; and 
constant field (Z-axis) modeling coefficient. 



Magnetic Field Modeling Equation 

Finally, the field modeling coefficients and field decay powers are used in the magnetic field modeling equation. 
This equation reconstructs the three components of the magnetic field on the X-Y plane of the magnetometers as a function 
of radial distance from the center of rotation (R) and the angle of rotation ( 6 ). 

Bj (R, 9) = DC j /r Aj + H,(C,.«) cos(nd)/R Bn ‘ +H,(^«)sin(n^)//?^ 

Assuming four triaxial magnetometers are used in the dipole moment measurements (i = 1 to 4), the components (for n=4) 
from the X-axis magnetometers are: 

CXj = (2A„/R? -6A 31 /Rf + 45A S1 /4R 7 -35A 71 /2Rf) 

SXj = (2B, , /Rf -6B 3) /Rf +45B SI /4R, 7 -35B 71 /2R?) 

The components from the Y-axis magnetometers are: 

CY, = - (B, j /Rf -3B 31 /2Rf +15B S1 /8R 7 -35B 71 /l6R?) 

SYj = (A„/R? -3A 31 /2Rf +15A 5 ,/8R 7 -35A 71 /l6R?) 

The constant component from the Z-axis magnetometers is: 

DCZj = -(A 10 /R? -3A 30 /2Rf +15A S0 /8R 7 -35A 70 /l6Rf) 

Where: 

An = X-axis dipole moment coefficient; 

A 3 i, A 5 i, A 7 i = Contributions to the fundamental component from the X-axis multi-pole (octapole, etc.) moments; 

Bit = Y-axis dipole moment coefficient’ 

B 3 i, B 5 i, B 7 i = Contributions to the fundamental component from the Y-axis multi-pole (octapole, etc.) moments; 

A io = Z-axis dipole moment coefficient; and 

A 3 o, A 50 , A 70 = Contributions to the constant component from the Z-axis multi-pole (octapole, etc.) moments. 

The Fourier analysis equations to measure the fundamental and constant components from rotational magnetic field 
data taken every 10° are: 

CX * = * f B x(tf,,0)cos<?c/0 - 

SXj f'XMsin BdO= 

CY , =i £B y (R„e)co*6d6 = tVl^Mcos* 

SY , =^[ X B y (R i ,0)sm0d0=^lZ B >M sin0 

dcz., =£ f b XR„e)de = j-XZ 

Where: 

B x ,B y ,B z = X-, Y- and Z-axis magnetic fields (nT); 

Ri = Magnetometer radial distance (m); 

0 = Rotation angle; and 

1 = Magnetometer number, 1 to 4. 

X Moment Near-Field Analysis 

The X Moment Fourier coefficient and Near Field matrices are as follows: 

CXj 2 (RjR,) 3 -6 (RJR,) 5 45 (RjRj -35 {RjR,) 9 

oCX, oCX, aCX, 4 aCX, 2 oCX, 

SY,' = (j RJR,) 3 -3 (RJR,) 5 1 5(RjR,y -35 {RjR,f 

oSY, aSY, 2 aSY, 8 aSY, 1 6aSY, 

Q Matrix column , column 2 column 3 column 4 



Solving for the Data Matrix results in the following: 

DM, =^Q Matrix x column, = Z(RJR,) 3 ^CXjiaCX,) 2 + SY l l(oSY l f) 

DM 2 =£q M atrix xcolumn, = £ (7?, /i?, ) 5 x (- 6CX, /(oCX , ) 2 -2SY i l2(dSY i f) 

DM, = YjQ Matrix x column^ = £ (R } /R,) 1 x {aSCX, /4(oCX, ) 2 + 1 557, /8 (aSY t ) 2 ) 

DM 4 = Y,Q Matrixxcolumn > = XW/^-) 9 x(- 35 C*i/ 2 (oC* ( ) 2 - 3 5S7,/ 1 6(o5^) 2 ) 

Solving for the Measurement Matrix results in the following: 

MM,, =^(coK) 2 =X(^iM) 6 x((2/oOf,) 2 +(l/oSJ') 2 ) 

MM 22 =XM« 2 ) 2 =XWM) 10 x((- 6 /oCA ' / ) 2 +(-3/2oSY i ) 2 ) 

MM 33 = 'Z(column 3 ) 2 =^(7?,/5,) 14 x ((45/4oCX,.) 2 + (15/8o57,) 2 ) 

MM44 =^](co/um« 4 ) 2 =^(i?,/5,) 18 x((-35/2oCX,) 2 + (-35/16a57,) 2 ) 

MM 12 = MM 21 = £ (co/wwin, ) x (co/uw« 2 ) = £ (i?, /tf, ) 8 x (- 12/ (oCX,) 2 -3/2(aSY,) 2 ) 

MM,3 = MM„ = 2 (column, ) x (column , )=£(*, //?, ) 10 x ^0/4(oCX i ) 2 + 1 5/8(o5nK ) 2 ) 

MM, 4 = MM 4I = £ (column , ) x (column , ) = £ (5, /i?, ) 12 x (- 70/2(oC7, ) 2 -35/16(oS^) 2 ) 
MM 23 = MM 32 =Y(column 2 )x (column , )=£(/?,//?, ) 12 x(- 270/ 4(oCJ7,) 2 -45/1 6(057 ) 2 ) 
MM 24 = MM 42 =^(co/w/n« 2 )x(co/wm« 4 ) =^(^/J?,.) 14 x(210/2(oCJ7,) 2 +105/32(oS7 j ) 2 ) 
MM 34 = MM 43 = ^ (column, ) x (column , )=£(/?,//?, ) 16 x (- 1 575/8(oC7, ) 2 - 525/1 28(o5^. ) 2 ) 

Thus, solving for the dipole moments requires the following matrix: 


DM, 


MM n 

mm 2I 

mm 3 , 

mm ai 


A, 

DM, 


MM n 

mm 22 

mm, 2 

MM n 

v 

A: 

DM, 


MM n 

mm 2 , 

MM„ 

MM t , 

A 

^51 

DM 4 


mm , 4 

mm 24 

mm 34 

mm 44 


A 7l 


For the first iteration dipole moment calculation (N=l), a lxl matrix is produced that reveals dipole fields only. Solving for 
the dipole moment: 

Moment= A,, xlOR 2 
Moment = (DMl/MM,,)x 10 R 2 

For the second iteration dipole moment calculation (N=2), a 2x2 matrix is produced that reveals dipole and octapole fields. 
Solving for the dipole moment: 


DM1*MM22-DM2*MM21 ,. n 

Moment = — — — - x 10 R 

MM1 1 * MM22 - MM12 * MM21 


3 

1 


For the third iteration dipole moment calculation (N=3), a 3x3 matrix is produced that reveals higher order multipole fields. 
Solving for the dipole moment: 


DM, x(MM 22 xMM 33 -MMj 3 xMM 32 )- 
MM 21 x (DM 2 x MM„ - DM 3 x MM 32 ) + 


Moment = 


MM 3 , x (DM 2 x MM 23 - DM, x MM 22 ) 
MM n x(MM 22 xMM 33 -MM 23 xMM 32 )- 


xlOR 2 


MMj, x (MM, 2 x MM 33 - MMj 3 x MM 32 ) + 


MM„ x(MM , 2 xMM 2! - MM , 3 xMM 22 ) 



Y Moment Near-Field Analysis 

Solving for the Y Moment Data Matrix results in the following: 


DM, ='£(RjR,y X^sxjiosxy -CYjicCY,) 1 ) 

DM 2 = £ (/?, /R, ) 5 x(- 6SX, /(pSX, ) 2 + 3CK, /2(oCY, ) 2 ) 
DM, = £ (RJR I ) 7 x(45&Y f /4(o5Y,) 2 -15C^/8(oC^) 2 ) 
DM 4 = x (-355Y,./2(o5Y,) 2 +35CF,/16(oCy,) 2 ) 


Solving for the Measurement Matrix results in the following: 

MM„ ='£(RjR l .) 6 x{(2/oSX i ) 2 +(l/aCY i ) 2 ) 

MM 22 = Y (R, /R, )'° x k~6/oSX i ) 2 + (-3/2 ocy, ) 2 ) 

MM„ =Z(^M) 14 x ((45/ 4o5Y,) 2 + (15/8oCr,) 2 ) 

MM,, xii-^/loSX,) 2 +(-35/l6oCY i ) 2 ) 

MM 12 = MM 21 = ^ (*, /*, )‘ x (- n/ioSXtf - 3/2(aCY f ) 2 ) 

MM„ =MM„ = £(/?, /tf^x (90/4(oSY,) 2 +l 5/8(00^ ) 2 ) 

MM, 4 =MM 4I =X(^M) 12 x(-70/2(oSY,) 2 -35/16(oCY,) 2 ) 
MM a = MM, 2 = £(/?,//?, ) 12 x(-270/4(a£¥,) 2 -45/16(oCl') 2 ) 
MM 24 = MM 42 = ^(RJR,)' 4 x(210/2(oSY,) 2 +]0S/32(oCY i ) 2 ) 
MM, 4 = MM 4 , = £(/?,//?,) 16 *(-1575/8(oSY,) 2 -525/128(oCT,) 2 ) 
The Y Moment is calculated using moment equations for the X-axis. 


Z Moment Near-Field Analysis 

Solving for the Z Moment Data Matrix results in the following: 
DM, = £ (R t /tf,) 3 x (- DCZ, /(ctDCZ, ) 2 ) 

DM, =^(/?,//?,) 5 x(3DCZ,/2(<tDCZ,) 2 ) 

DM, = x(-15DCZ,/8(<jDCZ,) 2 ) 

DM 4 =^(^,/^,) 9 x(35DCZ,./16(oDCZ,) 2 ) 

Solving for the Measurement Matrix results in the following: 
MM„=X(«,M) 6 x(1/oDCZ,) 2 
MM 22 =X(*>M) 10 x (-3/2ctDCZ ; ) 2 
M M„ = £(tf,/tf,) 14 x(15/8ccDCZ,) 2 
MM 44 =£(/?,//?,.)“ x (-35/1 6o£)CZ,) 2 
MM|j = MM 2 , = ^(i?,/*,) 8 x(- 3/2(000^ ) 2 ) 

MM„ = MM„ =X(^,M)'° x(l5/8(oDCZ,) 2 ) 

MM, 4 =MM 4 , =X(^ 1 /^) 12 x(-35/16(oDCZ i ) 2 ) 

MM 2 , = MM, 2 = £(/?,//?,) 12 x (- 45/1 6{aDCZ i ) 2 ) 

MM 24 = MM 42 = ^](/?,/^) 14 x(l05/32(o£»CZ j ) 2 ) 

MM,, =MM 4 , =^(/?,/i?,.) 16 x(-525/128(o£)CZ,.) 2 ) 

The Z Moment is calculated using moment equations for the X-axis. 


The oCX;, oSYj, oSXj, oCYj and oDCZj variables are moment constant term standard deviations derived from the 
average data values calculated based on five rotational runs. Experience has shown that only one rotational run is necessary. 
Since only one rotational run is used, the standard deviation values are equal to 1 and drop out of the dipole moment 
calculations used in the Lab VIEW rotational data analysis software. 

Discussion 

The magnetic field dipole moment modeling technique has been implemented on a microcomputer-controlled data 
acquisition and analysis system at the Goddard Space Flight Center Magnetic Test Facility. This system is described in [2] 
and consists of National Instruments (NI) hardware and NI Lab VIEW software. The hardware includes a NI plug-in data 
acquisition board inside the microcomputer (or a data acquisition card in a laptop) cabled to an external Signal Conditioning 
extensions for Instrumentation (SCXI) chassis with signal SCXI conditioning equipment. Lab VIEW is an object-oriented 
graphic programming code used for data acquisition, analysis and presentation. 

The basic assumption that Fourier coefficients have the same polarity for all magnetometers and decrease with 
increasing radius leads to certain spatial limitations of the technique. In the extreme near-field region, complex magnetic 
sources with higher order multi-pole moments can overshadow lower order moment contributions to the Fourier coefficients. 
This condition allows the possibility of polarity reversals with decreasing radius. Therefore, to maximize accuracy, the region 
to be modeled should be narrowly defined by the placement of in measurement magnetometers around complex magnetic 
sources. 


The possible calculation of total field, three-dimensional magnetic field models is offered by rotating the spacecraft 
and acquiring data for different values of the vertical angle <p, from 0 to 7t. This provides data for models on R-0 planes as a 
function of <p, allowing modeling on R-tp planes as a continuous function of 0. This would yield coefficients sufficient to 
provide modeling over all R, 0 and <p parameters; however, it is felt that most modeling requirements could be met with 
simple mapping on a single R-0 plane or over just a limited region of the vertical angle. 

Conclusions 

Test results obtained from modeling magnetic fields around several different spacecraft and flight experiments 
indicate that the magnetic field modeling technique, combined with a computer-controlled measure system, provides 
capability for complex magnetic dipole moment analysis. Modeling accuracy is sufficient to provide detailed analysis of 
complex magnetic sources, even where field levels are less than 1 .0 nT. 

Presentation of the magnetic dipole moment model can he achieved in real-time, using high-speed data acquisition 
and analysis augmented by a CRT and laser printer. Similar systems have been incorporated into existing magnetic 
measurement facilities and also could be designed as part of the spacecraft checkout platform on the Space Shuttle. Magnetic 
testing using the Shuttle would allow measurements to be made on spacecraft operating in true flight configurations, 
significantly increasing the reliability and utility of the data. 
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